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Neutrino oscillation experiments and direct bounds on absolute masses constrain neutrino mass 
differences to fall into the microwave energy range, for most of the allowed parameter space. As 
a consequence of these recent phenomenological advances, older constraints on radiative neutrino 
decays based on diffuse background radiations and assuming strongly hierarchical masses in the eV 
range are now outdated. We thus derive new bounds on the radiative neutrino lifetime using the 
high precision cosmic microwave background spectral data collected by the Far Infrared Absolute 
Spectrophotometer instrument on board of Cosmic Background Explorer. The lower bound on the 
lifetime is between a fewxlO^^ s and ~ 5 x 10^'' s, depending on the neutrino mass ordering and on 
the absolute mass scale. However, due to phase space limitations, the upper bound in terms of the 
effective magnetic moment mediating the decay is not better than ~ 10~* Bohr magnetons. We also 
comment about possible improvements of these limits, by means of recent diffuse infrared photon 
background data. We compare these bounds with pre-existing limits coming from laboratory or 
astrophysical arguments. We emphasize the complementarity of our results with others available in 
the literature. 
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I. INTRODUCTION 

The last decade has seen two great improve- 
ments in astroparticle physics: the wealth of in- 
formation on neutrino physics (for a recent re- 
view see e.g. and the impressive precision of 
new cosmological data (sec for example the latest 
WMAP team results 0]). The standard cosmologi- 
cal scenario predicts the existence of a diffuse back- 
ground of low-energy neutrinos, which has been of- 
ten investigated in the past to probe non-standard 
neutrino properties (for a few recent examples, 
see 0, 13,11, Ig). At this stage, it is meaningful 
to reassert the impact of non-standard physics in 
the neutrino sector on cosmological observables, or 
equivalently to re-cxplore the constraints provided 
by cosmology on exotic physics. 

In this paper we revisit the bounds on neutrino 
radiative lifetime coming from cosmology. Indeed, 
older constraints based on the diffuse Cosmic In- 
frared Background (CIB) and assuming strongly 
hierarchical masses in the eV range 0, Q (see 
also 0, [13, [HI for general reviews) are now out- 
dated and strictly speaking inapplicable. The neu- 
trino mass splittings squared provided by oscil- 
lation experiments and present upper bounds on 
the neutrino mass scale constrain neutrino mass 
differences to fall in the microwave energy range 
{E ^ 10~^ eV), in most of the allowed parame- 
ter space. Proper bounds must be derived using 
the high precision cosmic microwave background 



(CMB) data collected by the Far Infrared Abso- 
lute Spectrophotometer (FIRAS) instrument on 
board the Cosmic Background Explorer (COBE), 
which tested the blackbody nature of the spec- 
trum at better than 1 part in lO" [H, [H. The 
high precision of this measurement has been al- 
ready exploited to constrain some new physics sce- 
narios, producing deformations on the CMB spec- 
trum [l5| . In addition, a new estimate of the 
CIB flux (although less accurate than CMB) has 
been recentlyderived [ll| from the SPITZER tele- 
scope data [13]. Consequently, one can also infer 
updated bounds on the radiative lifetime when the 
unrcdshifted photon energy falls in this infrared 
range {E ^ 10~^ eV), i.e., in the limit of non- 
degenerate neutrino mass pattern. 

Here we perform such an analysis for radiative 
neutrino decays, and compare the new bounds ob- 
tained with limits coming from laboratory or astro- 
physical arguments. Although taken at face value 
these bounds are not competitive with the most 
stringent ones already available, we emphasize that 
they probe different combinations and/or regimes 
of the effective couplings describing the electro- 
magnetic properties of neutrinos, thus being com- 
plementary (rather then redundant) with respect 
to the others available in the literature. The plan 
of our work is as follows. In Sec. |TT]we summarize 
the relevant formalism, while in Sec. Illll we present 
the CMB data used and the bounds obtained. In 
Sec. II VI we briefly discuss the bounds coming from 
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the CIB flux. In Sec. |V]we comment our results 
and give the conclusion. 



II. RADIATIVE NEUTRINO DECAYS 

Let us denote by Vi the (active) neutrino fields 
respectively of masses m^, where i = 1,2,3. The 
radiative decay vt — > ^'j + 7 can be thought of as 
arising from an effective interaction Lagrangian of 
the form 

Ant = i^v'-aapi^itj + e^jj^y^ F°''^ + h.c. (1) 

where F"'^ is the electromagnetic field tensor, 
<^ap = [7a , 7/3] where 7^ are the Dirac-matrices and 
[. , .] is the commutator, Vi is the neutrino field of 
mass rrii, and and eij are the magnetic and 
electric transition moments usually expressed in 
units of the Bohr magneton jis- The convention 
to sum over repeated indices is used. In general. 
fiij and eij are functions of the transferred momen- 
tum squared q^, so that constraints obtained at a 
different arc independent. The radiative decay 
rate for a transition i —>■ j is written 
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In the following, we shall quote the bounds in 
terms of . We shall assume that the radiative 
decay rate is very low compared with the expansion 
rate of the universe; neither the cosmological evolu- 
tion or the primordial neutrino spectrum is affected 
by the additional coupling we are going to intro- 
duce. A posteriori, this is known to be an excellent 
approximation. For the same reason, we shall also 
neglect "multiple decays" (the daughter neutrino 
Vj constitutes a negligible fraction of the original 
Vi quasi-thermal population). We shall take our 
input data for neutrino mass eigenstate densities 
from the calculation performed in without any 
extra parameter, as non- vanishing chemical poten- 
tials. With present data, the latter are anyway 
constrained to be well below 0{1) [l^, so drop- 
ping this assumption would not change much our 
conclusions. 

From simple kinematical considerations it fol- 
lows that in a decay Ui i^j + ^ from a state of 
mass rrii into one of mass ruj < rui, the photon 
in the rest frame of the decaying neutrinos is thus 
monochromatic (two-body decay), with an energy 



(3) 



At present, the neutrino mass spectrum is 
constrained by the well-known values of the 
two squared mass splittings for the atmospheric 
{Amjj) and the solar (Am|) neutrino problems. 
We take their best-fit values and 2a ranges from 
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7.92(1 ±0.09) 

2.6(lt°:l^)xl0- 
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The remaining unknowns in the neutrino spectrum 
are the absolute mass scale (equivalently, the mass 
of the lightest eigenstate mi) and the mass hierar- 
chy. Namely, in normal hierarchy (NH) the mass 
pattern would be 



mi 



m2 



ma = \/ml + Am| + Amj^ ; (6) 
while in inverted hierarchy (IH) one would have 
mi , 



Am|^ ; 
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In the limiting case of normal hierarchy and 
mi = 0, the lightest neutrino for which a decay 
is possible has a mass m2 ~ 9 x 10"'^ eV and is 
thus non-relativistic for most of the universe life- 
time, namely in the redshift range z < 50. We can 
thus safely work in the approximation of all neu- 
trinos decaying effectively at rest. In this limit, we 
can also neglect the momentum distribution of the 
neutrino spectra. The formalism which would al- 
low one to generalize our results to the momentum- 
dependent case has been developed in [23|, which 
we address the interested reader for further details. 
However, the corrections arc small, of the order of 
powers of the neutrino temperature to mass ra- 
tios, and also vanishing in the limit of very long 
lifetimes. 

We shall discuss the limits on as a function 
of mi and for the two cases NH and IH. We shall 
vary the mass scale in < mi < 2 eV as allowed by 
the Mainz experiment on the ^H beta decay end- 
point [2l|. In this respect, we shall be conserva- 
tive: If neutrinos are Majorana particles the more 
stringent bound from Oz^/3/J searches apply, with an 
effective mass bound mpp < 0.81 eV Struc- 
ture formation, combined with other cosmological 
data, also constrains m^ 0]. Present cosmolog- 
ical bounds span the range ^ 0-2 ~ 2eV, 
0, US, m, Hi, m, ill depending on the data 
sets used and priors assumed. An upper limit of 
^■m^ ~0.6 eV (i.e. mi ~ 0.2 eV) is often con- 
sidered robust, and we shall report it for illus- 
trative purposes. Yet, one may circumvent the 



3 



OvfiP bound (e.g. with a Dirac neutrino) and 
significantly relax the most stringent cosmological 
bounds (for example with a conservative combina- 
tion of cosmological data sets and priors or with an 
exotic dark energy sector), so in the following we 
shall present our results up to the value mi = 2 cV. 

Let Fe be the present energy flux of photons 
with present energy E produced by neutrino decay. 
The differential energy flux ipE (energy flux Fe 
per unit energy and solid angle) is related to the 
differential number flux (the particle flux Fn 
per unit energy and solid angle) at present by 



d^FE 

dE dn 



= E 



dE dn 
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and it can be shown that, if the lifetime Ti of the 
neutrino of mass is much greater than the uni- 
verse lifetime it holds 1201 ^ 
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47r H(Z32) 4^ i?(23i) 47r H{z2i) 
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where rii ~ 113 cm ^ is the present number den- 
sity of the i— th neutrino in absence of decay, 
the Hubble function is (assuming, for simplicity, 
a flat cosmology) H{z) ^ i7o\/^^A/(l + z)^ + Ha, 
Hf) ~ 73 km s~-^ Mpc~^ being the present Hubble 
expansion rate, and CIm — 0.26 and CIa ~ 0.74 
respectively the matter and the cosmological con- 
stant energy density relative to the critical one. 
The dependence on energy enters implicitly via the 
quantities < = £ij / E — 1. 

In practice, to a very good approximation one 
can write a general equation of the kind 
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where, however, the meaning of the factors how- 
ever depends on the hierarchy. In NH, in the 
first two terms of the sum in Eq. ^ it holds 
Z32 cii Z31 = ZH, and one can identify zl ~ Z21, 

ri = r]i, = r;^^ -i- r;^2- ih, it is the 

last two terms of the sum in Eq. © which have 
Z31 ~ Z21 = zh, and using 712 ~ one can iden- 
tify Zl = Z32, Tl = rj^, rjj = r^^ + r],^. in 

both cases, we shall therefore express our bounds 
in terms of keeping in mind their slightly dif- 
ferent meaning for the two cases of NH and IH. 

In Fig. 1 we represent the unredshiftcd photon 
energy Sij from decaying neutrinos [Eq. (3)] as a 
function of the lightest neutrino mass eigenstate 



^ Note that Ti may be much shorter than the radiative 
Ufetime, which in most exotic models is dominated by 
invisible decays l27t |28|| . In the present work we are 
neglecting these cases. 
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FIG. 1; Unredshifted photon energy e from decaying 
neutrinos [Eq. (3)] as a function of the lightest neu- 
trino mass eigenstate mi, for the two neutrino mass 
splittings (L,H) in normal and inverted hierarchy. (See 
text for details) The horizontal band represents the 
energ y r ange of the CMB spectrum measured by FI- 
RAS hj. The GIB energy range is also shown. 



mi in the case of normal and inverted hierarchy, 
where the meaning of eL,H is clear from the pre- 
vious discussion. We also indicate by an horizon- 
tal band the energy range of the CMB spectrum 
(2.84 X 10"'' eV < E < 2.65 x IQ-^eV) measured 
by FIRAS [l2[. We also show the CIB range in the 
energy band above the FIRAS range up to (con- 
ventionally) 0.15 eV [H. For mi < 0.5 eV, the 
photon energy eh falls in the CIB range. 

For photons emitted at z = in the FIRAS 
range, the effect of radiative decays is most promi- 
nent and results in a feature on the CMB spectrum. 
Actually even if photons arc emitted at higher en- 
ergy the effect is still strong, since photons emitted 
at a redshift of a few enter the FIRAS spectrum be- 
cause of cosmological redshift; it is easy to check 
that one has thus some sensitivity to kh in the 
whole range for mi. However, as we will see, for 
mi < 0.1 eV a stronger (but less robust) limit can 
be obtained using directly the CIB data. 

On the other hand, for mi > 0.14 eV the pho- 
tons corresponding to the smaller splitting are 
falling in the radio band, below the frequency 
range probed by COBE, where measurements are 
more uncertain and thus one has no sensitivity to 
kl and the corresponding bound disappears. 



III. THE CMB BOUND 

To constrain the neutrino electromagnetic de- 
cay we use the COBE/FIRAS data for the exper- 
imentally measured CMB spectrum, corrected for 
foregrounds [l^. Note that the new calibration of 
FIRAS [l^l is within the old errors and would not 
change any of our conclusions. The = 43 data 
points ^^^^ at different energies Ei are obtained by 
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FIG. 2: Bounds on th and tl vs. mi, for the two cases 
of NH and IH. The regions below the sohd curves are 
excluded at 95 % C.L. The curves for th coincide in 
the two cases, although the definition of th is different 
(see text) . The dot-dashed line represents the limit on 
Th obtained from cosmic infrared background. 

summing the best-fit blackbody spectrum to the 
residuals reported in Ref. p^ . The experimental 
errors (Ji and the correlation indices pij between 
different energies are also available. In the pres- 
ence of neutrino decay, the original radiance (en- 
ergy flux per unit of solid angle) of the "theoretical 
blackbody" at temperature T 

<i>"iE,T)^^[cME/T)^iy' (11) 

vifould gain an additional term so that the intensity 
becomes 

= ^"{E,T) + ipEiKl^H,rni). (12) 
We then build the reduced chi-squared function 
1 ^ 

(13) 

where 

= $rP-$(K„T) (14) 
is the i-th residual, and 

<^t3 = PrjCr.CTj (15) 

is the covariance matrix. In principle, the parame- 
ter T entering initially in $"(£', T) needs not to be 
fixed at the standard value Tq = 2.725 ± 0.002 K 
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FIG. 3: Bounds on kh and kl vs. mi, for the two cases 
of NH and IH. The regions above the solid curves are 
excluded at 95 % C.L. The curves for kh coincide in 
the two cases, although the definition of hh is different 
(see text). The dot-dashed line represents the limit on 
kh obtained from cosmic infrared background. 

[Tsj . which is the best fit of the "distorted" spec- 
trum eventually observed now. The initial T before 
a significant fraction of neutrinos decays should 
thus be a free parameter, to be determined in the 
minimization procedure. Practically, however, the 
distortion introduced by the neutrino decay spec- 
trum is such highly non-thermal that a change in 
T can not be accomodated for any significant neu- 
trino lifetime: the constraints obtained fixing T are 
basically the same. 

Our results are reported in Fig. 2, where we 
represent the exclusion plot in the plane tl^h = 
(Tjf j^)^^ vs. mi, where the regions below the solid 
curves are excluded at 95 % C.L. For small values 
of mi the most stringent limit is > 4 x 10^" s 
in IH (slightly better than in NH case), while 
the bound on th is about an order of magnitude 
smaller, say r^j > 2 x 10^^ s, since for low rrii only 
photons produced by H decays at a redshift of 
few are in FIRAS range. On the contrary, for 
mi >0.14 eV, the bound on tl disappears, while 
the bound on th becomes more stringent, being 
th>5x 10^° s. Note that the "fuzzy" behaviour of 
the bounds is due to the sharp edge of the photon 
spectrum at E — Eh.l- when the photon energy 
embeds a new FIRAS bin, the function has a 
sharp discontinuity. 

In Fig. 3 we translate the plot of Fig. 2 in an 
exclusion plot in the plane mi vs. hl,h- Here 
the factor (Smfj/rrii)^ maps in a non trivial way 
the bounds in terms of kl,h- The regions above 
the solid curves are excluded at 95 % C.L. For 
the NH case, kl ^ 3 x 10~^ ^b, while in the IH 
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case, kl <3 x 10~^/iB. In agreement with our 
previous considerations, the bound on kl dis- 
appears for mi >0.14eV. On the contrary, the 
bound for kh is always present, and it corresponds 
to kh < 8 X 10 /iB apart for the degenerate re- 
gion, where it degrades down to 10~^ jiB or even 
more. Note also that typical cosmological upper 
bounds would already exclude the extreme degen- 
erate case. 



IV. THE CIB BOUND 

The Cosmic Infrared Background (CIB) is 
mainly the relic emission of the formation and 
evolution of the galaxies of all types at wave- 
lengths larger than a few microns. The spec- 
trum of the CIB is peaked around ~ 100 fim 
[E ~ 1.2 X 10~^ eV), thus just is in the energy 
range £h of photon from radiative v decays, for 
mi < 0.1 eV. Recently, a new estimate of the 
CIB flux has been established using the Spitzer 
Observatory data [l^. The measured CIB flux is 
^ciB ~ 24 nW m~^ sr~^. Using this number we 
can obtain a rough bound on th (and hence on 
k,h) simply requiring that the total energy flux of 
the photons coming from v decay does not exceed 
the CIB flux: 



i^E dE < $CIB 



(16) 



where we consider as lower limit of the CIB range 
the upper value of the FIRAS range, i.e. i?min = 
2.65 X 10^'^ eV. The bounds of th and kh obtained 
from Eq. p6p are shown respectively in Figs. 2 and 
3 by the dot-dashed line ^ . Although these bounds 
are stronger than those obtained by the FIRAS 
data in the same range of mi, we emphasize that 
they should be considered only as indicative. In 
fact, the CIB flux have still strong uncertainties 
and the precise spectral shape is essentially un- 
known (a factor ~ 3 of uncertainty should be ac- 
counted [H]). 

It is interesting to comment that, if it turns out 
that mi < 0.1 eV, an improvement on the bound 
on Th will clearly take advantage of a better mea- 
surement of the CIB flux and a more detailed 
knowledge of the astrophysical sources contribut- 
ing to it. Conversely, mi > O.leV would imply 
a signiflcant degree of neutrino clustering in large 
dark matter halos, the larger the mass the stronger 
the clustering [2^ [lOl . In turn, the expectation of 



^ For NH and mi < 10~^ eV also falls marginally in 
the CIB range. However, we have explicitly checked that 
the FIRAS constraint on is always stronger. 



overdensities would motivate analyses in the mi- 
crowave sky toward specific targets (like nearby 
galaxies or galaxy clusters), thus taking full ad- 
vantage of the spectral feature expected from neu- 
trino decay and looking for an angular-dependent 
enhancement over the CMB background. In spirit, 
this would be similar to what performed in the X- 
ray band when searching for signatures of sterile 
neutrino or axion decays (see e.g. [3ll, H^). Al- 
though a detailed treatment of these issues goes 
beyond the purpose of this work, it is worth not- 
ing that which one is the regime to consider will 
be basically answered by the KATRIN experiment 
on tritium beta decay [33|. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have revisited the bounds 
on the neutrino radiative lifetime coming from 
cosmology, deriving updated constraints from the 
high precision CMB spectrum data collected by 
the FIRAS instrument on board of COBE. We 
also compare these bounds with those obtained (by 
a back-of-the-envelope calculation) using the mea- 
surement of the flux of the Cosmic Infrared Back- 
ground, which although sometimes overrides those 
coming from the CMB data, should be considered 
only as qualitative. 

Previous cosmological bounds were either de- 
rived in a pre-COBE era Q or there were assumed 
higher masses and/or a strongly hierarchical mass 
spectrum, thus using the infrared background to 
derive the constraints Q. This has motivate us 
to re-evaluate the bounds within the presently al- 
lowed range of parameters suggested by neutrino 
oscillation physics and tritium endpoint experi- 
ments. Since it is customary to parameterize the 
neutrino electromagnetic decay via an effective op- 
erator of the kind reported in Eq. ((T|), it makes 
sense to translate the bounds (which actually are 
on the lifetime) into bounds on the parameters 
kh,l [see Eqs. (|2ll0p ]. These constraints are not 
better than k < 10^* /is which at first sight do 
not appear competitive with astrophysical limits 
neither with most of the laboratory bounds [s^ . 
Nevertheless, every experimental measure and ev- 
ery cosmological and astrophysical constraint has 
its own systematic uncertainties and its own rec- 
ognized or un-recognized loop-holes. In this sense 
it is certainly important to use many different ap- 
proaches to constrain fundamental neutrino prop- 
erties. 

In particular, the cosmological bound is based 
on the appearance of the daughter photons, and 
thus is very direct (modulo the underlying cos- 
mological assumptions). Therefore the bounds 
on th,l are completely independent from the un- 
derlying model that mediates the neutrino de- 
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cay. Each decay model must face with our direct 
bounds on T}j,l wich are the strongest attainable 
with direct livetime measurements. Moreover, it 
also probes the energy scale < 10"'^ eV^, unac- 
cessible to both laboratory experiments and stel- 
lar arguments. Other constraints have different 
features. Neutrino electromagnetic couplings are 
actually tightly constrained by energy-loss argu- 
ments in stars, in particular via the plasmon pro- 
cess 7* — > i/z/, to be Kij < 3 x 10^^^ fis fssl. Isoj. 
These are indirect limits, which strictly speaking 
apply to > keV^. Laboratory bounds are ob- 
tained via elastic v e scattering, where the scat- 
tered neutrino is not observed, and are at most at 
the level of 10~^° /is for the electron flavor (3^.[37j. 
The combinations of matrix elements and Hij 
that arc constrained by various experiments de- 
pend on the initial neutrino flavor and on its propa- 
gation between source and detector. Cancellations 
may occur in exotic cases [loj , and in any case the 
energy scale probed is q^ > MeV^. An additional 
motivation for independent checks is that models 
with a strong dependence of Kij{q'^) have been pro- 
posed [11]. Our bound seems to exclude extreme 
runnings of the effective couplings. Vice versa, if 
one assumes that the elements ctj and fitj are quasi 
energy independent, laboratory and astrophysical 
arguments exclude any possibility of phcnomcno- 



logically interesting radiative decays in cosmology; 
this is a priori surprising, given that cosmology in- 
volves the longest time intervals available and the 
lowest boosting factors for neutrinos. On the other 
hand, note that the cosmological bound may be 
easily violated by invoking a neutrino invisible de- 
cay much faster than the universe lifetime. It is 
intriguing to notice that, despite the impossibility 
to test this "nightmare case" in laboratory, cos- 
mology may probe to some extent those scenarios, 
potentially with important consequences for parti- 
cle physics as well [l^, HI] . 
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